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Using the SU(3) C lattice QCD with the Kogut-Susskind fermion at the quenched level, we study 
the quark-gluon mixed condensate g(qa IJA/ G IJA/ q) , which is another chiral order parameter. For each 
current quark mass, m q = 21,36 and 52 MeV, we generate 100 gauge configurations in the 16 4 
lattice with (3 = 6.0, and perform the measurement of g^o^G^q) at 16 points in each gauge 
configuration. Using the 1600 data for each m q , we find mo = g(qo^ v G^ v q) / (qq) — 2.5 GeV 2 at 
the lattice scale in the chiral limit. The large value of giqa^G^q) suggests its importance in the 
operator product expansion in QCD. 

PACS numbers: PACS: 12.38.Gc, 12.38.-t, 11.15.Ha 


I. INTRODUCTION 


The main signature of the non-perturbative nature of quantum chromodynamics (QCD) is its nontrivial vacuum 
structure, which is represented by various condensates, or vacuum expectation values. For instance, the quark 
condensate (qq) is a standard order parameter of spontaneous chiral symmetry breaking in QCD, and it determines 
properties of hadrons, especially the pion and the other pseudoscalar Nambu-Goldstone bosons. In the gluonic 
sector, the gluon condensate (G^G^") is an important quantity associated with the trace anomaly in QCD, and the 
topological susceptibility ( Q 2 ) is responsible for the large rf mass due to the U^(l) anomaly. Recently, the behavior 
of the various condensates at finite temperature/density is a subject of intensive research in the context of the QCD 
phase diagram, particularly the transition to the quark-gluon-plasma phase. 

Among various condensates, we emphasize here the importance of the quark-gluon mixed condensate g(qcr^G M „g) = 
9(qpq)■ First, the mixed condensate represents a direct correlation between quarks and gluons in the QCD 
vacuum. In this point, the mixed condensate differs from the above-mentioned condensates even at the qualitative 
level. Second, this mixed condensate is another chiral order parameter of the second lowest dimension and it flips the 
chirality of the quark as 

giqa^G^q) = g(q R (a^G^) 9l) + g(qL (o>G M!/ ) q R ). (1) 


Note here that the mixed condensate plays a relevant role in the operator product expansion (OPE) in QCD as the 
next-to-leading chiral variant operator. 

Also for the low-energy phenomena of hadrons, the mixed condensate is found to be important through the frame¬ 
work of the QCD sum rule djj, which connects the various condensates in OPE and the hadronic properties with 
the help of the dispersion relation. The condensates are determined phenomenologically so as to reproduce various 
hadronic properties systematically, considering the Borel stability of the sum rules 0. For instance, in the standard 
QCD sum rule, the nucleon mass m n and the delta mass toa are given in terms of (qq) and g(qa M „G ^q) as|l(4j] 
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where M denotes the Borel mass and Xn and Aa parameters in the QCD sum rule Hi Here, we have used the 
standard parameterization as 


mo = g{q^^ v G llv q)/{qq). (5) 

In the QCD sum rules, the value mg ~ 0.8 ± 0.2 GeV 2 has been proposed as a result of the phenomenological 
analyses iiHi Therefore, for the sum rule on A, the second OPE term, which is proportional to the mixed 
condensate, amounts to the same magnitude to the leading OPE term, if we take the Borel mass M equals to the 
typical baryon mass as lGeV. From these equations, one finds that the condensate giqcr^G^q) has large effects on 
the N -A splitting 0. 

The condensate g(q<J^ v G^ v q) is also important in the light-heavy meson systems |||, since the term run■ g{qcr^G^q) 
proportional to the heavy quark mass m# contributes significantly in OPE of the corresponding sum rule. Further¬ 
more, through the direct mixing of q, q and gluons, the mixed condensate g(qo tLV G directly contributes in the 

exotic meson systems m. ' "" 

Needless to say, it is desirable to estimate g(qa^G^q), not only by the phenomenological parameter fitting in QCD 
sum rules, but also by a direct calculation from QCD. For this purpose, the lattice QCD Monte Carlo simulation [ll| is 
a powerful tool. With this method, the condensates can be directly calculated from QCD, keeping the non-perturbative 
effect. However, in spite of the importance of giqcr^G^q), there was only one preliminary lattice QCD study for 
g{q^fuyG^ v q) done about 15 years ago 0 This pioneering study 0 gave an estimate TOq ~ l.lGeV 2 , but this result 
is not conclusive yet because the simulation was done with insufficient statistics using a small and coarse lattice: the 
authors used only 5 gauge configurations on the 8 4 lattice with (3 = 5.7, and calculated the condensates at only 1 
space-time point for each gauge configuration. 

Therefore, in this paper, we present the calculation of g(qcr ^G ^ v q) in lattice QCD with a large and fine lattice 
and with high statistics. We perform the measurement of glqa^G^q) as well as (qq) in the SU(3) C lattice at the 
quenched level. Since these condensates are chiral order parameters, we adopt the Kogut-Susskind (KS) fermion, 
which keeps the explicit chiral symmetry in the massless quark limit. We generate 100 gauge configurations and pick 
up 16 space-time points for each configuration to calculate the condensates. Therefore, we obtain 1600 data for each 
quark mass and each (3. We perform reliable estimate of the condensates with this high statistics. 

This paper is organized as follows. In Sec. on we explain our formalism to calculate the condensates. In Sec. EH we 
present the lattice QCD data, and discuss its reliability by performing several checks. Sec. II VI is devoted to summary 
and concluding remarks. 


II. FORMALISM 


In this section, we describe the formalism on the calculation of the condensates (qq) and g(q<J fJ , l/ G^q} in SU(3) C 
quenched lattice QCD. Note that, even without the dynamical quark effects, the quenched lattice QCD calculations 
have reproduced various hadronic properties in good agreement with empirical values 0 . Moreover, the character¬ 
istics of the quenched simulation are well under control owing to the accumulated knowledge. Therefore, it is worth 
performing the quenched lattice QCD calculation before proceeding to the full QCD calculation as a next step. 

The lattice QCD is formulated in terms of the link-variable [/^(s) = exp[— iagA^(s)\ on the lattice with spacing a, 
instead of the continuum gluon held A^(x). For the gauge sector, we adopt the standard Wilson action as 


^ = EE^ 1 

S fI>L' L 


—Re Tr U^(s) 


( 6 ) 


with (3 = 2NJg 2 and the plaquette operator U^ v (s) on the (//, ^)-plane, which is described by 

U^(s) = U^(s)U v (s + g)Ul{s + v)Ul{s). (7) 

For the fermion action, we adopt the KS-fermion. As the advantage of the KS-fermion, its action takes a simple form 
and preserves the explicit chiral symmetry in massless quark limit, m = 0. The latter property of the KS-fermion is 
desirable for our study, since both of the condensates (qq) and g(qa^ V G^q) are expected to be sensitive to explicit 
chiral symmetry breaking as chiral order parameters. 

We comment here on the other lattice fermions briefly. The domain-wall fermion would be attractive from the 
viewpoint of chiral symmetry. However, its simulation is much more expensive in comparison with the KS-fermion. In 
addition, there are ambiguities originating from the newly introduced simulation parameters such as the domain-wall 
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height. The Wilson and the clover fermions would not be appropriate for our purpose, because they have a serious 
disadvantage from the viewpoint of chiral symmetry. Specifically, the action for these fermions contains the term 

£e x sb oc q(s)[U^(s)q{s ++ Ul(s-/j,)q(s - fi)-2q(s)] , (8) 

which explicitly breaks chiral symmetry even for m = 0. Although this term vanishes in the continuum limit, the chiral 
order parameters inevitably suffer the nontrivial contamination from this unphysical term at finite lattice spacing. 
This uncontrollable contamination should be avoided. 

The action for the KS-fermion [l3} with the mass m is described by 

Sf = 2 [ u »( s )x{s +P) - Ul{s - n)x{s - p)] +maY l X(s)x(s), (9) 

S,fl S 

where \ an d X are Grassmann fields which have no spinor degrees of freedom, and is the staggered phase defined 
as i)fj, (s) = (-l) 51+ "’ +5 "- 1 , i.e., 

m(s) = 1, r/2 = (-l) Sl , % = (-l) Sl+S2 , 7?4 = (-l) Sl+S2+S3 . (10) 

In order to make the definition of the sign of g(qcr tlv G tlu q) unambiguous, we note here that the definition of the 
continuum covariant derivative is D^ — igA corresponding to the definition of = e~ la9AlJ -. 

In this formalism, the quark field q is introduced as an SU(IV/ = 4) spinor field. The explicit relation between the 
quark field q and the spinless Grassmann field \ is understood in the following way. When the gauge field is set to 
be zero, the quark field q is expressed by 

4 0*0 = ( T phf X(x + P), ( 11 ) 

P 

T P = 7r7 2 P2 7 3 P3 7f, (12) 

where p = (pi, p 2 , P 3 , p±) with p^ £ {0,1} runs over the 16 sites in the 2 4 lrypercube. The indices i and / denote the 
spinor and the flavor indices, respectively. When the gluon field is turned on, we insert additional link-variables in 
Eq. ITT!) in order to respect the gauge covariance. 

The evaluation of the condensates amounts to the following expressions as 

a 3 {qq) = [ sf ( x ’ x )\ > ( 13 ) 

/ 

a 5 g(qv^ G ^q) = i sf ( x ^ x ) a ^ G ^] > ( 14 ) 

/ AC 1 ' 

where the trace “Tr” is taken with respect to both the spinor and the color indices, and x ) denotes the Euclidean 

quark propagator for /-th flavor as 

S f (y,x) = (q f (y)q f (x)). (15) 

In terms of x an( l X-Selds, the flavor-averaged quark condensate is rewritten on the lattice as 

a 3 {qq) = Tr [ r ^ r p + p)x{x + p))] ■ (16) 

p 

The corresponding diagram is shown in figure ^ 

On the other hand, the flavor-averaged quark-gluon mixed condensate is given by 

a 5 g(q<J^G^q) = Tr t U ±n,±»( x + P) r p' r p (x( x + p')x{ x + p)) cr^ G^(x + p) ] , (17) 

P,*' p 

p' = p ± p ± v, 

where the sign ± is taken such that the sink point (x + p') = {x + p ± p ± v) belongs to the same hypercube of the 
source point (x + p). Here, in order to respect the gauge covariance, we have used in Eq. ra 

(*^) = ^ (p^) U±i/ {x dl /i) U±i/ (x) {x zh Z/) j , 


( 18 ) 
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FIG. 1: The diagrammatic representation of (qq) in terms of the spinless Grassmann field \ on the lattice. The solid curve 
with the arrow denotes the propagation of x- 



FIG. 2: The diagrammatic representation of the two ingredients in g(qo ^vG ^q) in terms of the spinless Grassmann field x 
and the gluon field on the lattice, (a) The left diagram shows the propagation of x> with the insertion of gauge links. The 
solid curve with the arrow denotes the propagation of x, and the curly lines denote the inserted gauge link U M ,„. (b) The right 
diagram shows the gluon held strength Gjf*, where each loop of the curly line denotes a plaquette operator. 



where we use the definition of U- I _ l (x ) = U^(x — fj,). 

On the gluon held strength G^ v , we adopt the clover-type definition on the lattice, 

G» = ^ E AA Tr i xA i U ^ s ) + U »~^ s ) + U-^{s) + £/■_„(*)} - \ A W - u}\ , (19) 

A 

where \ A (A = 1, 2, • • • , 8) denotes the color SU(3) Gell-Mann matrix normalized as Tr(A‘ 4 A B ) = 2 5 AB . In hgure|3 
we show the diagrams corresponding to Eqs. and m - 
In the continuum limit, Eq. m leads to 


GjSS(a) 


+ 0(a 2 ) 


( 20 ) 


It is worth mentioning that this dehnition has no 0(a) discretization error. On the other hand, in Ref.|ll2l|. the authors 
adopted a simple insertion of the gluon held strength, 


1 


a b g(qa^G^q) = - ^ q f {s)a^U t _ l (s)U v {s + g)q f (s + g + v), 


( 21 ) 


which contains O(a) error. Although both of the dehnitions, Eqs. il 71) and CD, coincide to the giqcr^G^q) in the 
continuum limit, our dehnition of G will give less systematic errors in the actual lattice simulations with finite 
lattice spacing a. 


III. THE LATTICE QCD RESULTS 
A. Lattice QCD results for g{qo llv G llu q) 

We calculate the condensates (qq) and g^o^vG^q) using the SU(3) C lattice QCD at the quenched level. The Monte 
Carlo simulation is performed with the standard Wilson action for (3 = 5.7, 5.8 and 6.0 on the 8 4 ,12 4 and 16 4 lattice, 
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TABLE I: The lattice QCD parameter (3 = 2 N c /g 2 and the lattice size used in the simulation. The lattice spacing a, the 
physical volume V and the adopted values of the current quark mass ma are also listed for each f3. As for the current quark 
mass, the corresponding physical values are m = 21,36 and 52 MeV from the left. 


(3 

lattice size a [fm] 

V[fm 4 ] 


ma 


5.7 

8 4 

0.19 

(1-5) 4 

0.0200 

0.0350 

0.0500 

5.8 

12 4 

0.14 

(1.7) 4 

0.0147 

0.0258 

0.0368 

6.0 

16 4 

0.10 

(1-6) 4 

0.0105 

0.0184 

0.0263 


TABLE II: The numerical results of a 3 {qq) and a 5 g{qa llv G llv q) for various ma in SU(3) C lattice QCD with (3= 6.0 and 16 4 . 
The last column denotes their values in the chiral limit obtained by the linear chiral extrapolation. 

ma = 0.0263 ma = 0.0184 ma = 0.0105 chiral limit 
a^Jqq) -0.04240(16) -0.03247(15) -0.02212(16) -0.00872(17) 

a 3 g{qa^G^q) -0.01882(15) -0.01498(14) -0.01088(14) -0.00565(14) 


respectively. The pseudo-heat-bath algorithm is adopted for the update of the gauge configuration. After 1000 sweeps 
for the thermalization, we pick up 100 gauge configurations for every 500 sweeps. The lattice unit a is determined 
so as to reproduce the string tension a = 0.89GeV/fm 0 . In Table m we summarize the lattice parameters for the 
gauge configuration. We note that the physical volume V is roughly the same for the three calculations with different 

(3. 

We use the quark mass parameter, m = 21,36 and 52 MeV, which correspond to ma = 0.0105, 0.0184 and 0.0263 
for /3 = 6.0, respectively. Also for f3 = 5.7 and 5.8, we use the same values of the physical quark mass m. The 
corresponding values of ma are also tabulated in Table |I] 

In the determination of the Euclidean propagator {x(yx)x( x x)) ■ we s °l ve the matrix inverse equations iteratively 
using the CG, BiCGSTAB and MR algorithms, until the residual error r 2 becomes small enough to satisfy r 2 < 
10 -8 (/3 = 6.0) or r 2 < 10 -1 °(/3 = 5.7, 5.8). By checking the differences of the results among these algorithms, 
we confirm the numerical errors are smaller than the statistical errors for all /3. For the Grassmann %-field, the 
anti-periodic condition is imposed. The dependence on the boundary condition will be discussed later. 

In the KS-fermion formalism, the source point x x = x + p of the y-field is taken to be on the hypercubic site around 
the physical source point x. We take 16 physical space-time source points x in each gauge configuration as follows: on 
the lattices with the volume (2n) 4 = 8 4 ,12 4 and 16 4 , we take x = (x\, X 2 , £ 3 , X 4 ) with x M G {0, n} in the lattice unit. 
For each physical space-time point x , we take the sum over p in the hypercube, corresponding to the flavor and spinor 
contractions. For each j3 and m, we calculate the flavor-averaged condensates according to Eqs. lP|l and and 
average them over the 16 physical space-time points and 100 gauge configurations. Statistical errors are calculated 
using the jackknife error estimate. 

Figure 0 shows the values of the bare condensates a 3 (qq) and a 5 g(qa l _ lv G^q) against the quark mass ma. We 
emphasize that the jackknife errors are almost negligible, due to the high statistics of 1600 data for each quark mass. 
From figure 0 both (qq) and g{qcr^G^q) show a clear linear behavior against the quark mass to. This feature is 
also found for (3 = 5.7 and 5.8. Therefore, we fit the data with a linear function and determine the condensates in the 
chiral limit. The results are summarized in Tables mi mu and hyi 


TABLE III: The numerical results of a 3 {qq) and a 5 g(qa l _ l ^G fll/ q) for various ma in SU(3) C lattice QCD with (3= 5.8 and 12 4 . 
The last column denotes their values in the chiral limit obtained by the linear chiral extrapolation. 

ma = 0.0368 ma = 0.0258 ma — 0.0147 chiral limit 
a 3 {qq ) -0.08031(35) -0.06667(35) -0.05159(36) -0.03271(38) 

a 5 g{qa^G^q) -0.03883(32) -0.03375(32) -0.02783(31) -0.02064(32) 
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FIG. 3: The bare condensates a 3 {qq) and a 3 g(qa IJU G i,, v q) plotted against the quark mass ma at /3 = 6.0. The dashed lines 
denote the best linear extrapolations, and the cross symbols correspond to the values in the chiral limit. The jackknife errors 
are hidden in the dots. 




TABLE IV: The numerical results of a 3 (qq) and a 5 g(qa llv G l _i, v q) for various ma in SU(3) C lattice QCD with (3= 5.7 and 8 4 
The last column denotes their values in the chiral limit obtained by the linear chiral extrapolation. 


ma = 0.050 ma = 0.035 ma = 0.020 chiral limit 


a°{qq ) -0.12346(68) -0.10788(73) -0.09017(80) -0.06833(90) 

a 5 g{qa^G^q) -0.06200(53) -0.05695(56) -0.05064(60) -0.04327(66) 
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TABLE V: The lattice results of a 3 {qq) and a 5 g(qa l _ l ^G IJ , v q) with (3= 6.0 and 16 4 in the case of the periodic boundary condition 
for the fermion field. This calculation is done for the check of the boundary effect and the finite volume artifact. 

ma = 0.0263 ma = 0.0184 ma = 0.0105 chiral limit 
a 3 (qq) -0.04236(16) -0.03240(16) -0.02200(16) -0.00854(17) 
a 3 g(qa^G^q) -0.01880(15) -0.01493(14) -0.01078(14) -0.00551(14) 


B. Check on the systematic uncertainty 

In this section, we check the reliabilit y of our lattice QCD results. We first consider the finite volume artifact. As 
indicated by the Banks-Casher formula Jlfj , 


!im lim (qq) = -tt——, (22) 

m—>0 V —too V 

for the spectral density p( A) of the Dirac operator, the total volume V should be large enough before the quark mass 
goes to zero. In order to estimate this finite volume artifact, we carry out the same calculation imposing the periodic 
boundary condition on the Grassmann fields x and X-, instead of the anti-periodic boundary condition, keeping the 
other parameters same. If the lattice total volume is too small, the quark propagates over the total volume, and the 
results would be sensitive to the boundary conditions. Thus, the difference will indicate ambiguity from the finite 
volume artifact. 

We show in Tabic IVI the lattice results with the periodic boundary condition for (3 = 6.0. Comparing with Table ITT1 
one finds that the difference is only 1% level. The similar results are obtained for (3 = 5.7 and 5.8. Therefore, we 
conclude that the physical volume V ~ (1.6 frn) 4 in our simulations is large enough to avoid the finite volume artifact. 

We next consider the discretization error. As an advanced feature of the KS-fermion, the discretization error begins 
from 0(a 2 ) on the lattice spacing a. This is because O(a) errors cancel with each other when the average over the 
SU(4) flavor is taken. Therefore, there is no 0(a) error originating from the quark propagator in Eqs. 11311 and US- 
On the other hand, there is an ambiguity coming from a particular choice of the gauge link in Eq. ra, which is 
introduced to respect the gauge covariance. This ambiguity can be checked by changing the definition of g(Q<J t ivG liI/ q), 
adopting a different path which connects the source point (x + p ) and the sink point ( x + p') in Eq. 11711 . Specifically, 
instead oiU^. v in Eq. <d, we examine the other product as 

U±n,±v{x) = U± ll {x)U± v {x ± p), ( 23 ) 


and thus define giqcr^G^q) by 


giqa^G^q) = Tr U±n,±v( x + P) r p' r p M x + p')x(x + p)) a^ G l ^{x + p) 


(24) 


P 


where p' = p±p±v and the sign ± is taken as before. We perform the calculation for (3 = 5.7 and ma = 0.050. In this 
case, the lattice spacing a is largest in our simulations, and therefore the discretization error is expected to be larger 
than the other cases for (3 = 5.8 and 6.0. At each gauge configuration, we check the difference of the results between 
the choice of Eq. m and Eq. (|23ll. Since the typical difference is about 1%, we conclude that the discretization error 
is small enough, which confirms the reliability of our lattice results. 


C. Determination of m 2 = g{qa^G^q) / (qq) 

The values of the condensates in the continuum limit are to be obtained after the renormalization. To this end, 
the lattice perturbation theory has often been used, although it is afflicted with uncertainty originating from the non- 
perturbative effect. In principle, the non-perturbative renormalization scheme is desirable, which, however, requires 
a lot of computational power [lfij . Therefore we seek for another way which can reduce this uncertainty. Here we 
provide the ratio = g{qa^G^q) / {qq) , where some of the uncertainties are canceled with each other. In particular, 
this ratio is free from the uncertainty from the wave function renormalization of the quark. As a consequence, the 
results become more reliable with less uncertainties. In addition, the dependence of mg on the lattice spacing is 
weakened to a 2 , while (qq) and g(qa ^„G ^q) are proportional to a 3 and a 5 , respectively. We note that m 2 itself has 
the following physical meaning. In the QCD sum rule, g(qcr^G^q) generally appears as the chiral variant term next 










to (qq) in OPE. Therefore, m g is usually used without referring to the absolute value of g (qcr ^G ^ v q) itself, and thus 
it represents the level of importance of g(qa M „g) in OPE. 

Now, we present the result of the ratio ml using the bare results in SU(3) 0 lattice QCD. We adopt the results at 
(3 = 6.0, since its lattice spacing is the finest in our calculations. We find in the chiral limit 

m 2 0 = g(qo^G^q)/(qq) ~ 2.5 GeV 2 (/3 = 6.0), (25) 

at the lattice cutoff scale as a -1 ~ 2 GeV. This large value of Wq suggests the significance of the mixed condensate 
in OPE. Although we do not include renormalization effect, this result itself is determined very precisely. 


IV. SUMMARY AND DISCUSSIONS 

We have studied the quark-gluon mixed condensate g(q<T liV G tiV q) using SU(3) C lattice QCD with the Kogut-Susskind 
fermion at the quenched level. First, we have emphasized that the mixed condensate is one of the key quantities in 
various quark hadron physics, especially in the baryon sector such as the N -A splitting. In spite of its importance, the 
lattice QCD studies of this quantity have been limited to only one preliminary study for 15 years. Recently, due to the 
progress in lattice QCD Monte Carlo calculations, it becomes possible to calculate this mixed condensate with much 
better statistics on a finer and larger lattice. For each quark mass of m q = 21,36 and 52 MeV, we have generated 
100 gauge configurations on the 16 4 ,12 4 and 8 4 lattice with (3 = 6.0, 5.8 and 5.7, respectively. We have performed 
the measurements of g(q<T IJlI/ G tll/ q) as well as (qq) at 16 physical space-time points in each gauge configuration. Using 
the 1600 data for each m q , we have found = g(qa^G^q) / (qq) ~ 2.5 GeV 2 in the chiral limit at the lattice 
scale a -1 ~ 2 GeV corresponding to (3 = 6.0. We have checked that the systematic and statistical errors are almost 
negligible. Therefore, the value of ml at the lattice scale has been well determined in this calculation. 

Finally, we compare our result with the standard value employed in the QCD sum rule. To this end, we change 
the renormalization point from /x ~ 7r/a to /x ~ 1 GeV corresponding to the QCD sum rule. Following Ref. 0 , we 
first take the lattice results of the condensates as the starting point of the flow, then rescale the condensates using 
the anomalous dimensions evaluated in a perturbative manner. We have the following rescaled condensates as 0, 


(qq) 

_ f Oi s (g) 
\a s (ir/a ), 

^ -4/bo 

<W>L /a> 


(26) 

g(qcr^G^q) 

( <x a (n) 

2/(3b 0 ) 

g(qcr^G^q) 


(27) 

tj. \a s (ir/a), 

) 

tv/ a 


where we use the one-loop formula, a s (/x) = bp in^/A^ —) Aqcd = 200 —300MeV and b 0 = (11/3)AT C — (2/3 )Nf. 

(The anomalous dimension given in Refs. 00 for g(qa^q) is different from Ref. 0 . However, this difference 
does not change our semi-quantitative analysis here.) For the case of the quenched lattice QCD, we adopt Nf = 0. By 
using our bare lattice QCD results at (3 = 6.0, we obtain m 2 1 = g^cr^G^q)/(qq)\^ ~ 3.5 — 3.7 GeV 2 at /x = lGeV, 
from (qq)\ ~ -(0.0477 - 0.0506) GeV 3 = -(0.36 - 0.37GeV) 3 and g(qa fi vG fJ . v q)\ ~ -(0.176 - 0.177) GeV 5 . 

Comparing with the standard value of ml = 0.8 ± 0.2 GeV 2 in the QCD sum rule, our calculation results in a rather 
large value. Note that the instanton model has made a slightly larger estimate as ml ~ 1.4 GeV 2 at /x ~ 0.6 GeV [2Cj . 
For a more definite determination of TOq, the renormalization procedure should be performed more carefully, which is 
also expected to improve the value of (qq) simultaneously. In principle, the non-perturbative renormalization scheme 
is most desirable, which would, however, require a significant calculation cost |16|. 

We again emphasize that the mixed condensate g(qa^G^q) plays very important roles in various contexts in quark 
hadron physics. Hence, it is preferable to perform further studies. In particular, the dynamical quark effects would 
be nontrivial, since the mixed condensate includes quark field. The thermal effects are also interesting in relation 
to chiral restoration, because the mixed condensate is another chiral order parameter. Actually, we are in progress 
with these two studies on the lattice 0|. Considering the RHIC project, it becomes more and more important to 
understand the nature of finite temperature QCD. Therefore, it is quite desirable to determine thermal effects on the 
condensates with lattice QCD in understanding the finite temperature QCD. 
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